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ABSTRACT 

WL16 is a unique member of the embedded young stellar population in the nearby p Ophiuchi 
cloud core: its extended, high surface brightness disk is visible only at mid-infrared wavelengths. 
We present diffraction-limited images, from 7.9 to 24.5 pm, of WL16 acquired at the Keck II 
telescope. We take advantage of the ~0'/3 angular resolution of the mid-infrared images to derive 
physical parameters for the central object by self-consistently combining them with available near- 
infrared spectroscopy, point-spread-function fit photometry, and pre-main-sequence evolutionary 
tracks. We find the central star to be a 250 L 0 , 4 M 0 , Herbig Ae star, seen through foreground 
material of the p Oph cloud core that provides an extinction of Ay = 31 ± 1 magnitudes. 
WL16’s disk is detected through all nine observed passbands, not only those four which sample 
PAH emission features. We confirm, therefore, that the emitting particles are composed of both 
polycyclic aromatic hydrocarbons (PAHs) and very small (5-100 A) graphitic grains (VSGs). 
The disk size as observed through the four PAH filters is 7"x3"5 , corresponding to a disk 
diameter of ~ 900 AU. The disk’s major axis is at a position angle of 60° ± 2° and viewed at an 
inclination angle of 62.2° ± 0.4° to our line-of-sight. Our derived inclination angle is in excellent 
agreement with the previously inferred inclination for the inner disk (R < 30 R 0 ) from kinematic 
modeling of the near-infrared spectral lines of CO. We can distinguish structure within the PAH 
disk at unprecedented resolution. We confirm a resolved (1"5 diameter) core component at 7.9 
and 8.8 pm, due to emission from positively charged PAHs. An enhancement in the emission 
at 12.5 pm at the disk’s edges is found for the first time, and signals the presence of larger 
(> 50-80 carbon atoms) PAHs and/or more hydrogenated PAHs than those found in the bulk of 
the disk. We find a disk asymmetry, observed at all nine mid-infrared wavelengths, at projected 
radii l"-2'.'5 (corresponding to 125 AU < r < 300 AU) from the central source. 
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1. Introduction 

WL16 was first discovered via near-infrared pho¬ 
tometry (Wilking & Lada 1983) to be a member 
of the deeply embedded cluster currently forming in 
the highest extinction core of the nearby p Ophiuchi 
cloud (Barsony et al. 1997, and references therein); 
it is undetectable at optical wavelengths. WL16 was 
identified as a few xlO 5 yr old, late-stage proto- 
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star based on its rising infrared (2-60 /am) spectral 
slope, and was assumed to consist of a remnant in¬ 
falling dust and gas envelope surrounding a central 
protostar-(-disk system (Lada 1987; Adams, Lada, 
and Shu 1987; Wilking, Lada, and Young 1989). 

Subsequent observations conflicted with the clas¬ 
sification of WL16 as a late-stage protostar, how¬ 
ever. WL16 lacks, by two orders of magnitude, the 
amount of millimeter flux expected from a low-mass 
young stellar object (YSO) at such an early evolu¬ 
tionary stage, and shows no evidence of a molecular 
outflow (Cabrit & Andre 1991; Andre & Montmerle 
1994). WL16 is the only source in Ophiuchus known 
to exhibit the entire suite of polycyclic aromatic hy¬ 
drocarbon (PAH) emission features (Tanaka et al. 
1990; Hanner et al. 1992; Deutsch et al. 1995; Natta 
& Kriigel 1995). It has become clear that PAH emis¬ 
sion accounts for the appearance of WL16’s mid- 
infrared spectral features rather than what had pre¬ 
viously been interpreted as a deep silicate absorption 
feature (Hanner et al. 1992), such as is generally as¬ 
sociated with the presence of a cool circumstellar 
envelope in late-stage, low-mass (< 2 M 0 ) proto¬ 
stars. 

WL16 is also the first YSO in which the pres¬ 
ence of a hot (1000K < T < 6000K), dense, neu¬ 
tral, gaseous, Keplerian disk was inferred via de¬ 
tailed modeling of its 2.3 pm CO overtone emission 
line shapes (Dent & Geballe 1991; Carr et al. 1993; 
Chandler et al. 1993, 1995). The inner disk (3— 
30 Rq) around WL16 is modeled as having a high 
inclination angle (> 48°) with a fairly low accretion 
rate from the disk to the star of M ~ 2 x 10” 7 Mq 
yr _1 (Najita et al. 1996a). Velocity-resolved spec¬ 
troscopy of the Bi '7 hydrogen line toward WL16 
provides independent evidence for the low accretion 
rate infalling gas at spatial scales smaller than that 
of the inner disk radius (Najita, Carr, and Tokunaga 
1996b; Hartmann, Hewett, and Calvet 1994). 

What makes WL16 truly unique, however, is 
the large extent (8" along its major axis) of its 
clliptically-shaped mid-infrared emission, as re¬ 
vealed by previous imaging studies (Deutsch et 
al. 1995; Emerson et al. 1996; Moore et al. 1998). 
Spatially-resolved mid-infrared spectroscopy of this 
elongated dust structure has confirmed the presence 
of both very small grains (VSG’s) and polycyclic 
aromatic hydrocarbons (or PAH’s) throughout (De¬ 
Vito & Hayward 1998). 

The presence of mid-infrared emission out to such 
large extents from a central source is not expected 
from “classical” or “standard equilibrium-emitting” 
silicate and amorphous carbon or “MRN” grains 


(Mathis, Rumpl, and Nordsieck 1977). These “stan¬ 
dard” grains range in size from 10 nm to > 250 nm, 
have a size distribution, n(a) oc a” 3 ' 5 , and are gener¬ 
ally assumed to be large enough to come into equilib¬ 
rium with the surrounding radiation field (although 
see Lynch & Mazuk 2000 for a thorough discussion of 
the shortcomings attendant upon this assumption). 

By contrast, the so-called graphitic “Very Small 
Grains” (VSGs), which range in size from 0.4 nm 
< a < 10 nm with a size distribution n(a) oc a” 4 
(Draine & Anderson 1985) and the PAHs, which 
contain 25-40 carbon atoms, or PAH clusters, 
containing 50-500 carbon atoms (Siebenmorgen 
& Kriigel 1992), are heated and emit in a non¬ 
equilibrium fashion. They are therefore capable 
of emitting in the mid-infrared, at large distances 
from a central illuminating source at (seemingly) 
“high temperatures” not in equilibrium with the 
surrounding radiation field. 

Detailed radiative transfer modeling of the spec¬ 
tral energy distribution of physically realistic, spher¬ 
ically symmetric dust shells surrounding Herbig Ae 
stars shows that whereas the mid-infrared (~ 2- 
20 /.mi) emission from “classical” grains is negli¬ 
gible, the continuum emission from VSGs in this 
wavelength range accounts for all of the observed 
mid-infrared excesses (above photospheric) for these 
sources (Natta, Prusti, and Kriigel 1993). These 
authors also note that VSGs do not emit in the J 
band (1.25 /im), and only very little in the H band 
(1.65 /im) for the same models in which VSGs con¬ 
tribute essentially all of the observed mid-infrared 
continuum emission. Therefore, the slope of the 
spectrum in the 1.25-1.65 /im continuum is effec¬ 
tively determined by the underlying stellar radiation 
field and of the absorption and scattering properties 
of the “standard” grains, even in the presence of 
VSGs. 

Whereas VSGs essentially determine the mid- 
infrared continuum emission in the circumstellar en¬ 
vironments of Herbig Ae/Be stars (Herbig 1994), 
the effect of adding PAH emitters results in emis¬ 
sion into several, well-defined solid-state features, 
with minimal change in the underlying shape of the 
mid-infrared continuum, as compared with mod¬ 
els computed in the absence of PAHs. When a 
UV/visible photon from the central source inter¬ 
acts with a PAH lattice, the photon energy is redis¬ 
tributed into stretching and bending modes of C—H 
and C=C bonds; these modes re-emit energy in 
broad, solid-state features. These well-known PAH 
features are centered at 3.29 /am (the v = 1 —> 0 
C—H stretch), at 6.2 and 7.7 /am (C=C stretch 
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modes), 8.6 pm (C—H in-plane bend), 11.3 pm 
(C—H out-of-plane bend), and 12.5 /tm (C—H bond 
resonance) (Tokunaga et al. 1991; Siebenmorgen & 
Kriigel 1992; Brooke, Tokunaga, and Strom 1993; 
Hanner, Brooke, and Tokunaga 1995). The region 
between 10 and 20 pm is dominated by a resonance 
between 10.5 pm < A < 13.5 pm (encompassing 
the 11.3 and 12.5 pm features) due to C—H out- 
of-plane bending. Although the approximate wave¬ 
length of one resonance is about 11.3 /im, due to a 
mono-C—H vibration, this wavelength can actually 
be longer, depending on hydrogenation (the num¬ 
ber of H atoms attached to the edge carbon rings). 
Some continuum emission is also produced by the 
overlapping wings of these PAH features. 

The excitation of mid-infrared emission by PAHs 
and VSGs out to such large physical distances from 
the central object as is observed in WL16, requires 
both a fairly luminous source of optical/UV pho¬ 
tons and optically thin paths from the central ob¬ 
ject to the mid-IR emitting regions. The unique¬ 
ness of the mid-IR emitting structure observed in 
WL16 amongst all of the known embedded sources 
in p Oph, combined with the fact that many inter¬ 
mediate mass (2 M 0 < M <8 M 0 ) YSO’s show 
PAH-feature emission (Waters & Waelkens 1998), 
leads one to re-consider the central source’s proper¬ 
ties. 

In order to shed further light on both the na¬ 
ture of the central object and its surrounding cir- 
cumstellar material, we have obtained images of 
WL16 at nine distinct wavelengths—through filters 
which variously either include or exclude PAH fea¬ 
ture emission and span the 8-25 pm atmospheric 
windows—all at unprecedentedly high spatial res¬ 
olution (~ O'/S). In addition, we have made use 
of previously published near-infrared photometry, 
recently published spectroscopy, and up-to-date 
pre-main-sequence tracks for intermediate stellar 
masses, to re-examine the nature of the central ob¬ 
ject. 

2. Observations 

We observed WL16 with JPL’s mid-infrared cam¬ 
era, MIRLIN (Ressler et al. 1994), at the visitor port 
of the Keck II telescope on UT 14 March 1998 and 
27 January 1999. The sky was clear and dry on both 
nights (t 225GHz ~ 0.06 and ~ 0.04, respectively). 
MIRLIN employs a Boeing HF-16, 128x128 pixel, 
Si:As impurity band conductor detector array, and 
produces a plate scale of 0'/138 per pixel (17"5 field- 



Fig. 1. - MIRLIN filter transmission curves along 
with a continuum-subtracted spectrum of WL16 
(DeVito & Hayward 1998). The filter curves are 
reproduced from the manufacturer’s data taken at 
77 K. The filters at 7.9, 9.7, and 11.7 pm have had 
their baselines shifted by one for clarity. The 9.7 
and 10.3 pm filters sample the “continuum” reason¬ 
ably cleanly, while the 7.9, 8.8, 11.7, and 12.5 pm 
filters each sample a PAH feature. 


of-view) with the f/40 chopping secondary mirror at 
Keck II. 

Background subtraction was performed by chop¬ 
ping the telescope secondary mirror 8" in a north- 
south direction, then nodding the entire telescope 
30" east-west, completely off the source, in order 
to remove residual differences. Observations were 
performed through the six 10 pm “silicate” filters 
(7.9, 8.8, 9.7, 10.3, 11.7, and 12.5 pm) and the 
17.9 /im filter on 14 March 1998, and through the 
20.8 /im and 24.5 pm filters on 27 January 1999. 
All of these filters have a passband of AA/A « 10%, 
except for 24.5 pm where AA/A ~ 3%. Plots of 
the 9 filter passbands are shown in Figure 1, along 
with the continuum-subtracted mid-infrared spec¬ 
trum of WL16 (DeVito & Hayward 1998), to indi¬ 
cate how the PAH emission features were sampled. 
Note that the well-known, broad silicate features, 
due to the Si-0 stretching and bending modes, at 
9.7 and 18 pm, respectively, were also sampled. 

Total on-source integration times through each of 
the “silicate” filters were approximately 30 seconds— 
200 coadded chop pairs of roughly 80 msec duration 
in each of the two beams. Total on-source inte¬ 
gration times through the longer wavelength filters 
were 2.0, 2.7, and 2.4 minutes through the 17.9, 
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Fig. 2.— Images of WL16 at all observed mid-infrared wavelengths. North is up, east is left for all images. 
The pixel scale is 0.138 arcsec/pix, giving the frames a roughly 13" x 9" field of view. The PSF calibrators 
were /3 Leo at A < 13 /jm, and a Hya at longer wavelengths, and are shown in the insets at each wavelength. 
Contours in each panel represent the flux density in Jy arcsec , and are spaced by one magnitude. The 
lowest contour level in all images starts at 50 mJy arcsec -2 ; subsequent contour levels are at 0.125, 0.315, 
0.792, 1.9, and 5.0 Jy arcsec -2 . The FWHMs of the PSFs at each wavelength are the diffraction-limited 
values of 0717 at 7.9 /tm, 0719 at 8.8 /mi, 0721 at 9.7 /tm, 0722 at 10.3 /tm, 0'/25 at 11.7 /tm, 0'/27 at 12.5 /tm, 
0739 at 17.9 /tm, 0745 at 20.8 /tm, and 0753 at 24.5 /tm. The “intensity profile marker” in the 12.5 /tm image 
shows the axis along which the brightness profiles of Figure 3 are plotted. 


20.8, and 24.5 /tm filters, respectively. The primary 
photometric standard and point spread function 
calibrator was the A3V star, (3 Leo, which has a 
magnitude ranging from 1.91 to 1.84 from 7-13 /tm. 
The calibrator at 17.9 and 20.8 /tm was a Hya, a 
K3III star with a magnitude of —1.49; the standard 
at 24.5 /tm was (3 Lib, a 2.84 magnitude B8V star. 
Consistency checks were performed with a CMa 
(mag = —1.39) and cr Sco (mag ~ 2.40); the latter 
proved to be an easily-resolved binary with 0745 
separation. Atmospheric extinction was corrected 
by observing the calibrators at several different air- 
masses. The resulting coefficients were found to be 
quite low (< 0.1-0.2 magnitudes per airmass). 


3. Results 

3.1. Mid-Infrared Images 

The mid-infrared images of WL16 are presented 
in Figure 2. It is clear that the circumstellar struc¬ 
ture surrounding WL16 appears brighter with re¬ 
spect to the central source in the filters that encom¬ 
pass the PAH features, relative to its appearance in 
the filters that exclude PAH feature emission over 
the 7 13 /mr region. Nevertheless, the mid-infrared 
emission is still well-resolved even in the filters that 
exclude PAH feature emission. 

The extent of the emission derived from the im¬ 
ages that include the PAH features is 7"x375 at 1% 
of the peak level. This projected source size corre- 
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Fig. 3.— Intensity profiles along the disk’s major axis at all observed wavelengths. In each panel, the 
solid curve is the intensity profile at the indicated wavelength, the dotted curve is the PSF profile at that 
wavelength, and the dashed curve is the scaled profile at 11.7 /tm (normalized at a point in the “shoulder” 
1.5 arcsec northeast from the center). 


sponds to 880 x 440 AU for an assumed distance 
of d = 125 pc (de Geus, de Zeeuw, and Lub 1989; 
de Geus 1992; Knude & Hog 1998). The position 
angle of the major axis, measured east from north, 
is 60° ± 2°. 

Intensity profiles along the major axis of WL16 
(indicated on the 12.5 /tm image of Figure 2) are 
presented for each wavelength in Figure 3. The 
central, bright core is resolved at the shortest (7.9 
and 8.8 /an) wavelengths, but is unresolved at wave¬ 
lengths longer than this. We also find that the con¬ 
trast of the central peak to the disk emission varies 
with wavelength (see Col. 5 of Table 1), with the 
highest PSF/disk contrast occurring in the 9.7 and 
10.3 /mi images and the lowest PSF/disk contrast 
in the 24.5 /im image. 

3.2. Mid-Infrared Color Ratios 

In Figure 4, we present three mid-infrared color 
ratio images of the WL16 disk: the 7.9/8.8 /im flux 
ratios in the top left panel, the 8.8/11.7 /an flux 
ratios in the middle left panel, and the 12.5/11.7 /an 
flux ratios in the bottom left panel. We find no 
significant color gradient in the entire 7.9/8.8 /an 
flux ratio image. This is in contrast to the other 
two flux ratio images presented in Figure 4. The 
8.8/11.7 n m image (and in the 7.9/11.7 /im image 
not shown here) shows significantly enhanced values 


within the inner ~ 1"5 diameter “core” (see also 
the corresponding panels of Figure 3). This inner 
core flux ratio enhancement is more quantitatively 
illustrated in the corresponding cross-cut panel of 
Figure 4, which shows the 8.8/11.7 /tm flux ratio at a 
peak value of ~ 1.75 within the inner core, declining 
to a relatively constant value of ~ 0.95 outside this 
core. By contrast, the 12.5/11.7 /tm flux ratio image 
shows little variation in the inner (< 2") regions; 
but shows a small enhancement towards the disk’s 
edges. This flux ratio rises from a value of w 0.85 in 
the inner regions to values > 1.1 at the disk’s outer 
edges. 

3.3. Photometry 

We present new mid-infrared and far-infrared 
photometry for WL16 in Table 1. To compare our 
photometry with previous ground-based results ( e.g. 
Lada & Wilking 1984; Moore et al. 1998), we have 
performed false aperture photometry with an 8" di¬ 
ameter software aperture. In general, the agreement 
between our 8-arcsec aperture photometry and that 
of Lada & Wilking is quite good in the mid-infrared. 
The random errors in our flux measurements are 
about 3% at 8-13 /tm, 7% at 17.9 and 20.8 /tm, 
and 30% at 24.5 /tm. Our fluxes agree with those 
previously published to 2 ct at 8.8 /tm, 3er at 10.3 /tm, 
and < la at 11.7 and 12.5 /tm. We find some dis- 
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Fig. 4.— Color ratio images of WL16 and corresponding cross-cuts taken along the disk’s major axis 
for selected wavelength pairs: Top panels: 7.9/8.8 fim, Middle panels: 8.8/11.7 /im, and Bottom panels: 
12.5/11.7 /im. North is up, east is left for each image. The pixel scale is 0.138 arcsec/pix, giving a 13"x9" 
field of view. The images were divided on a pixel-by-pixel basis and are presented unsmoothed. On the right 
hand side, we plot the corresponding flux ratios along the major axis as a function of angular distance from 
the central source. The ratio images were rotated by 30°, then the median and standard deviation were 
computed from a 7 pixel (1.0 arcsec) wide strip centered on the major axis and plotted accordingly. 
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crepancy, however, between our 9.7 /jm flux (2.2 Jy) 
and the previously published 9.7 /im flux (1.6 Jy) of 
Lada & Wilking, though this may not be significant 
due to the strong, telluric ozone absorption feature 
contained within this band. 

We also performed point spread function fit 
(PSF-fit) photometry of the central peak in each 
image of WL16 to better determine the flux from 
just the central source, since the spatial resolution 
of the mid-infrared images allows a clean separa¬ 
tion of the central peak from the surrounding neb¬ 
ulosity. To accomplish this, a PSF image derived 
from observed standard stars at each wavelength 
was shifted, scaled, and offset to fit the central 
source in WL16. This same method was applied 
to near-infrared images of WL16 to derive PSF-fit 
near-infrared photometry, as well (Ressler 1992). At 
all wavelengths, the PSF-fit flux contains just pho- 
tospheric light plus any warm dust within the PSF 
(< 40 AU). Values of these PSF-fit fluxes are also 
tabulated in Table 1. 

We present the spectral energy distributions 
(SEDs) for photometry measured through different- 
size apertures for WL16 in Figure 5. The data 
have been separated into PSF-fit points (triangles 
connected by a solid curve), 8 arcsec aperture pho¬ 
tometry points (squares and stars connected by a 
dotted curve), and the 45 arcsec resolution HIRES- 
processed IRAS fluxes (connected by dot-dashed 
lines). In the large (8-arcsec) aperture data, the 
steep drop in flux approaching the wavelength of 
9.7 /im from either side, was at one time mis¬ 
interpreted as a very deep silicate absorption fea¬ 
ture. The depth of this silicate absorption feature 
is significantly reduced in the curve connecting just 
the PSF-fit data points, underlining the important 
contribution of PAH feature emission to the larger 
aperture photometry. 

In the near-infrared, the flux values found within 
an 8-arcsec aperture (Wilking & Lada 1983), are 
only slightly higher than the flux values derived 
through PSF-fitting, in contrast to the large ob¬ 
served difference between the PSF-8" aperture pho¬ 
tometry values at all mid-infrared wavelengths, both 
through filters that include and those that exclude 
PAH features (see Figure 5). This behavior is as ex¬ 
pected of VSG emission: providing only a marginal 
contribution to near-infrared continuum emission, 
but being the dominant source of mid-infrared con¬ 
tinuum emission. This result also justifies the use 
of our near-infrared PSF-fit fluxes to derive the ex¬ 
tinction and intrinsic central source luminosity (see 
§4.1 below). 



Fig. 5.— The observed spectral energy distribu¬ 
tions (SEDs) of WL16. Triangles connected by a 
solid curve represent measurements derived from 
PSF fitting (filled triangles are near-infrared data 
from Ressler 1992; open triangles are mid-infrared 
data from this work). Fluxes derived from our 
aperture photometry (using 8-arcsec software aper¬ 
tures) are shown by open squares connected with a 
dashed curve. Data from Wilking & Lada (1983) 
and Lada & Wilking (1984) are plotted with as¬ 
terisks connected by a dotted curve. (The curve 
is broken between 4.8 and 8.8 /im since they ob¬ 
tained no data at 7.9 /im.) Crosses connected by the 
dot-dashed lines represent newly determined IRAS 
HIRES fluxes (this work). The dramatically increas¬ 
ing flux levels with increasing aperture observed at 
all mid-infrared wavelengths, including bandpasses 
that exclude PAH feature emission, attest to the 
presence of both spatially extended PAHs and VSGs 
in the circumstellar environment of WL16. 

4. Discussion 

4.1. Nature of the Central Source 

The spectral type of the central source of WL16 
has recently been determined from high signal- 
to-noise, moderate resolution, near-infrared spec¬ 
troscopy to lie in the range B8 to A7, corresponding 
to a photospheric effective temperature of 8,000 K 
< T e ff < 12,000 K (Luhrnan & Rieke 1999). We 
combine this constraint on the central object’s ef¬ 
fective temperature with other observational data 
and theoretical pre-main-sequence tracks in a self- 
consistent manner to re-determine both the esti¬ 
mated extinction to the central source and its intrin¬ 
sic bolometric luminosity, mass, and radius. The low 
millimeter flux (Andre & Montmerle 1994; Motte 
et al. 1998), lack of near-infrared excess (Moore 
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Table 1 

Photometry of WL16 


Wavelength 

H 

PSF Fitting 
Flux 

(Jy) 

Aperture Phot. 
Flux a 

(Jy) 

Prev. Publ. 
Flux 

(Jy) 

Disk/PSF Ratio 

(Aper.-PSF)/PSF 

Refs. 

1.2 

0.0037 


0.0040 


(1,2) 

1.65 

0.059 


0.071 


(1,2) 

2.2 

0.45 


0.48 


(1,2) 

3.4 

1.17 


1.36 


(1,2) 

3.8 

1.22 




(1) 

4.8 

2.18 


2.02 


(1,3) 

7.9 

1.49 

19.13 


11.8 


oo 

00 

0.65 

7.89 

8.3 

11.1 

(3) 

9.7 

0.35 

2.18 

1.6 

5.2 

(3) 

10.3 

0.30 

1.97 

1.8 

5.5 

(3) 

11.7 

0.49 

7.07 

7.1 

13.4 

(3) 

12.5 

0.57 

6.30 

6.25 

10.1 

(3) 

17.9 

0.28 

4.53 


15.2 


20.8 

0.25 

3.19 

4.8 b 

11.8 

(3) 

24.5 

0.19 

5.67 


28.8 


12 


16.0 (HIRES) 

18.8 


(4) 

25 


14.4 (HIRES) 

<40 


(4) 

60 


129 (HIRES) 

202 


(4) 

100 


223 (HIRES) 



(4) 

1300 



<0.006 


(5) 


a All quoted fluxes are for an 8" diameter aperture, except for HIRES fluxes which were 
determined from a 45" diameter aperture after local background subtraction. 

b Broadbancl Q measurement- AA/A ~ 40% 


References. —(1) Ressler (1992), (2) Wilking & Lada (1983), (3) Lada & Wilking (1984), 
(4) Young, Lada, & Wilking (1986), (5) Motte, Andre, & Neri (1998) 












et al. 1998), and the fact that UV/optical radia¬ 
tion is able to illuminate such a large (900 AU) 
VSG/PAH emitting region in the immediate vicin¬ 
ity of WL16, all lead to the conclusion that the 
circumstellar material surrounding WL16 is op¬ 
tically thin at these wavelengths. Furthermore, 
computations of the continuum spectral energy dis¬ 
tributions of flared and geometrically thin accre¬ 
tion disks have shown that high accretion rates 
(M > 10 -6 M 0 yr -1 ) are required to produce op¬ 
tically thick inner disks around HAeBe stars (Hart¬ 
mann, Kenyon, and Calvet 1993). For comparison, 
the upper limit for the inner disk accretion rate for 
WL16 is < 2 x 10 -7 M 0 yr _1 (Najita et al. 1996a). 

If the circumstellar environment of WL16 is opti¬ 
cally thin in the near-infrared, then the observed ex¬ 
tinction along the line-of-sight towards WL16 must 
arise in the intervening p Oph cloud, which acts as 
a foreground screen. We note that this is a radically 
different physical source picture than one in which 
all of the energetic UV/visible source photons are 
assumed to be intercepted by the immediately sur¬ 
rounding circumstellar material, and whose energy 
is subsequently re-radiated by this same circumstel¬ 
lar material at longer, infrared wavelengths. 

In order to determine the value of the intervening 
extinction, we fix the photospheric temperature of 
the central source at 9000 K. We include only the 
near-infrared photometry in the fit, since the near- 
infrared fluxes are essentially uncontaminated by 
any “classical” grain or VSG emission from circum¬ 
stellar material. We then fit a reddened blackbody 
(modeling the extinction with the Mathis (1990) 
curve for molecular clouds) with only two free pa¬ 
rameters, Ay and stellar radius, to determine the 
extinction through the cloud. 

The resulting extinction, Ay = 31, falls within 
the range of previously derived values. A strict up¬ 
per limit can be derived from the line-of-sight beam- 
averaged (55" FWHM) 13 CO observations, from 
which the extinction through the entire cloud in the 
direction of WL16 is estimated at Ay = 70 mag with 
a 50% error (Wilking, priv. comm.). At the other 
extreme, allowing for a wide range of intrinsic NIR 
source colors, a possible range of 2.0 < Ak < 2.9 
was found by Carr et al. (1993). A value for Ay = 31 
towards WL16 was previously published, derived 
from its H—K color by Wilking & Lada (1983). The 
most recent published determination for the extinc¬ 
tion towards WL16 is an Ay = 37, determined using 
the NIR photometry of Greene & Young (1992), and 
assuming an intrinsic A star photosphere (Moore et 
al. 1998). 


We argue that our determination of Ay = 31 ± 1 
for the extinction towards WL16 is the most strin¬ 
gent to date. Figure 6 shows data and a model 
that have all been dereddened for an intervening 
Ay = 31 screen using the Mathis (1990) extinc¬ 
tion curve. The solid curve plotted in Figure 6 
shows the resulting spectral energy distribution for 
a T = 9000 K blackbody photosphere with a ra¬ 
dius of 6.5 R 0 at d = 125 pc. The triangles repre¬ 
sent the dereddened PSF-fit data points; the aster¬ 
isks represent the corrected near-infrared data from 
Wilking & Lada (1983). Note that the residual 
silicate absorption feature apparent from the solid 
curve joining the mid-infrared, PSF-fit photometry 
in Figure 5 completely disappears from the corre¬ 
sponding extinction-corrected PSF-fit data plotted 
in Figure 6. This is especially striking in view of 
the fact that the mid-infrared fluxes were not in¬ 
cluded in the fit used to determine the extinction to 
the source. The dereddened, PSF-fit, mid-infrared 
data points do lie somewhat above the plotted pho¬ 
tospheric model curve, however. This slight, remain¬ 
ing mid-infrared “excess” may signal the presence of 
VSGs within < 40 AU of the central source (within 
the PSF). Nevertheless, the goodness of the fit sug¬ 
gests that the residual “dip” apparent in the mid- 
infrared, PSF-fit data of Figure 5 is caused by the 
intervening cloud’s silicate absorption, and not by 
PAH emission. 

Strong constraints on the possible values of the 
derived extinction can be made via the two require¬ 
ments of 1) matching the shape and magnitude of 
the dereddened near-infrared flux data to an appro¬ 
priate blackbody, and 2) eliminating the residual 
9.7 pm. silicate absorption feature from the curve 
joining the PSF-fit fluxes in Figure 5. If Ay were 
greater than ~ 32, the dereddened data would show 
a silicate emission feature in the mid-infrared and 
the near-infrared points would lie on a curve steeper 
than a blackbody. For Ay < 30, fits consistent with 
a 9,000 K blackbody emitter would become problem¬ 
atic. In fact, the derived Ay is relatively insensitive 
to the assumed temperature: variation from 6,000 
to 12,000 K produces only a 1 magnitude change in 
Ay since the infrared data are well on the Rayleigh 
portion of the blackbody curve. Therefore, we can 
confidently assert that the range of foreground ex¬ 
tinction falls within Ay = 31 ± 1. 

We can now arrive at a new determination of the 
bolometric luminosity of WL16, simply by integrat¬ 
ing under the photospheric curve plotted in Figure 6. 
We find a bolometric luminosity of 250 L 0 for an 
assumed distance of 125 pc. This is much higher 
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Fig. 6. - The solid curve represents the de-reddened 
spectral energy distribution (SED) of a T e // = 
9000 K, L* = 250 L 0 , i?* = 6.5 R 0 blackbody 
(photosphere) that represents the central illuminat¬ 
ing source of WL16, seen through an extinction of 
Ay = 31 ± 1 (see text). Corresponding de-reddened 
PSF-fit fluxes are indicated by filled triangles (near- 
infrared) and open triangles (mid-infrared); de- 
reddened fluxes measured through 8-arcsec aper¬ 
tures are indicated by stars (near-infrared). 

than previously published values for the luminos¬ 
ity of WL16, which range from 10 L 0 to 22 L 0 
for d = 160 pc (Young, Lada, and Wilking (1986); 
Wilking, Lada, and Young (1989)). This discrep¬ 
ancy is resolved when one takes into account the 
very different assumptions that went into the dif¬ 
ferent luminosity determinations. In the absence of 
evidence to the contrary, previous workers had as¬ 
sumed that all the extinction towards WL16 was lo¬ 
cal and circumstellar, and that all of the source pho¬ 
tons were absorbed and re-radiated locally at mid- 
and far-infrared wavelengths. However, in view of 
subsequent multi-wavelength photometry and spec¬ 
troscopy, the new source picture that emerges is one 
in which the extinction is provided by a foreground 
screen, remote from the source, so that the absorbed 
photons are reradiated well outside our (or anyone 
else’s) beam (see Figure 7 for the schematic cartoon 
of the source geometry). 

Our new bolometric luminosity determination 
of 250 L 0 , combined with the spectroscopically 
determined effective temperature determination of 
9,000 K, allows placement of WL16 on the H-R di¬ 
agram for the first time. Comparison of the posi¬ 
tion of WL16 on this diagram with the pre-main- 
sequence tracks of Palla & Stahler (1999) yield 


a source mass of ~ 4 M 0 and a source age of 
1.9 x 10 5 Yr. The Siess et al. (2000) PMS tracks 
also give a source mass of ~ 4 M 0 , but an age 
of 1.2 x 10 6 Yr. Previous determinations for the 
central source’s mass in WL16 were made before 
pre-main-sequence tracks for more massive objects 
(> 2.5 M 0 ) were available, coupled with an erro¬ 
neous underestimate (10-20 L 0 ) of the central ob¬ 
ject’s true luminosity. 

4.2. A Mid-IR Emitting Disk in WL16 

Previous mid-infrared imaging of WL16 shows 
an extended emission region elongated along an ap¬ 
proximately NE-SW axis (Deutsch et al. 1995; De¬ 
Vito & Hayward 1998; Moore et al. 1998). This 
emission has variously been interpreted as originat¬ 
ing from a nebula, a torus, a bipolar cavity, or a disk. 
We consider a “nebula” (meaning an H II region) to 
be unlikely given the lack of narrow spectral fea¬ 
tures in the near- or mid-infrared, and the fact that 
the central star has a much later spectral type than 
B3 (Luhrnan & Rieke 1999). It is also not a reflec¬ 
tion nebula given that scattering is inefficient in the 
mid-infrared, and the density necessary to show any 
scattering would require an optically thick disk; we 
have already argued that WL16’s emission is opti¬ 
cally thin. A torus is unlikely since the intensity 
profile does not indicate a “hole” in the inner re¬ 
gions of the disk, and a torus would be difficult to 
maintain physically. A bipolar cavity can be ruled 
out since the outflow to create such a cavity was 
searched for, but not found (Cabrit & Andre 1991). 
We propose that the most plausible interpretation of 
our multi-wavelength imaging dataset favors the no¬ 
tion that the mid-infrared emission originates from 
a geometrically thin, equatorial disk. 

Under this assumption for the source geome¬ 
try, the minor-to-major axis ratio is found to be 
0.466±0.007, measured from the images as observed 
in the four filters that include the PAH emission fea¬ 
tures. This ratio corresponds to an inclination angle 
of 62.2° ± 0.4° of the disk to our line-of-sight, if we 
assume the disk to be intrinsically circular and to 
appear elongated only because of projection effects. 
This value agrees well with the inclination angle 
(i = 60°) deduced for the inner disk (R < 30 R 0 ) 
of WL16 from kinematic modeling of the hot gas 
(T = 5000 K) emitting in the near-infrared CO lines 
(Dent & Geballe 1991; Carr et al. 1993; Chandler et 
al. 1993, 1995). 
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Ionized PAHs 



Fig. 7. A schematic showing the proposed configuration of WL16. The entire WL16 system (star plus 
large disk) is viewed through a portion of the p Oph cloud, which acts as an Ay — 31 foreground screen. The 
dust disk around the star is effectively optically thin, so that UV radiation propagates through the entire 
disk. The inner portions of the disk experience a UV field sufficiently strong to ionize a substantial portion 
of the PAHs, while a decrease in the field allows a higher fraction of large PAHs to exist at the periphery of 
the disk. 


4.3. Properties of the Mid-IR Emitting Disk 

It has already been suggested that the mid-IR 
emitting dust surrounding WL16 consists of at 
least two components, aromatic hydrocarbons and 
VSGs (Moore et al. 1998). The present observa¬ 
tions strongly support this suggestion. The disk of 
WL16 is easily visible through all the mid-infrared 
filters, not just those which contain PAH emission 
features, although the disk does emit most strongly 
through those filters which encompass PAHs (Fig¬ 
ure 2). Similar results were found by Moore et al. 
(1998), upon comparison of their images of WL16 
through filters which contain PAH emission features 
at 8.2 and 11.3 pm with their 10.2 pm continuum im¬ 
age, which is PAH-feature-free. Our results extend 
this finding to a greater number of wavelengths. 

Simple modeling of PAH and VSG emission in the 
mid-infrared from spherical distributions of optically 
thin dust surrounding an 8,000 K (12,000 K) cen¬ 
tral source, suggests that 1% (3.5%) of the starlight 
would be absorbed by PAHs and 15% (20%) by 
VSGs. However, these values for the PAHs are 
highly dependent on the assumed absorption cross- 
sections, and could be as high as 16.6% (17%) for an 
8,000 K (12,000 K) central object (Natta & Kriigel 


1995). For WL16, we find that ~ 1% of the 9,000 K 
central source’s luminosity, or 2.6 L 0 , are emitted 
by the PAHs and VSGs. At the very least, this 
low value, relative to the models, rules out spher¬ 
ical distributions of PAHs and VSGs around the 
source. For comparison, previous authors have esti¬ 
mated ~ 0.6 Lq to be contributed by the 7.7-13 pm 
PAH band emission within a 5" beam (Hanner et al. 
1992), scaled to d= 125 pc, and 1.2 L 0 , integrated 
over the entire disk extent (Moore et al. 1998). 

4-3.1. Very Small Grain (VSG) Constituent 

Figure 2 shows the extended appearance of WL16 
at 9.7, 10.3, 17.9, 20.8, and 24.5 pm; passbands 
which exclude PAH feature emission. Since neither 
PAHs nor the “standard” MRN grains can account 
for this emission, as we argue below, the continuum 
emitter responsible for this extended mid-infrared 
emission in the disk of WL16 must be in the form 
of the so-called very small grains (VSG’s) (Natta & 
Kriigel 1995; DeVito & Hayward 1998; Moore et al. 
1998). It has been shown that such VSGs surround¬ 
ing HAeBe stars are excited by a much broader spec¬ 
trum of photons than the PAHs and can reprocess 
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stellar radiation into the mid-IR (Natta, Prusti, and 
Kriigel 1993). 

We have already pointed out (see §4.1) that the 
systematic mid-infrared excess flux that remains af¬ 
ter de-reddening the PSF-fit photometry requires 
the presence of VSG’s on scales < 40 AU (see Fig¬ 
ure 6). Similarly, emission from VSGs is required 
to explain the large spatial extent of WL16’s disk 
at continuum wavelengths (R ~ 450 AU). This is 
quantitatively illustrated by Figure 8, where inten¬ 
sity cross-cuts of a model disk, consisting of regular 
ISM grains, are compared with the actual, observed 
intensity profiles of WL16’s disk. The model disk in¬ 
tensity profiles were kindly provided by Dr. Barbara 
Whitney, using a 3-D Monte Carlo code described 
in Wood et al. (2002), and includes both scattering 
and emission by standard interstellar grains (Kim et 
al. 1994) at mid-infrared wavelengths. The model 
shown is for a 0.01 Mg flared disk (the maximum 
disk mass allowed by millimeter observations), in¬ 
clined at 60° to our line-of-sight, at a distance of 
125 pc, illuminated by a central source of 250 L 0 
and T e ff = 9000 K. The disk density structure is 
as described in Wood et al. (2002), except the in¬ 
ner radius is set by the dust destruction radius for a 
T eff = 9000 K central source, and the disk’s outer 
radius is determined by the observations presented 
here. Note that emission from standard, interstellar 
grains at large radii from the central source would 
be undetectable, according to the models, in stark 
contrast to the observed emission. Standard, inter¬ 
stellar grains cannot, therefore, be responsible for 
the appearance of WL16’s disk at large radii through 
the filters which exclude PAH emission. 

Additional qualitative proof of this assertion is 
provided in Figure 3. A comparison of the 9.7, 
10.3, 17.9, 20.8, and 24.5 pm flux profiles (none of 
which should have significant PAH emission) with 
the fiducial 11.7 pm profile suggests a rather con¬ 
stant ratio as a function of distance from the cen¬ 
ter, despite any differences in the absolute surface 
brightness values. This also holds true in the fully 
two-dimensional case, as well, where ratios of any 
pair of images, including the 10.3 pm/20.8 pm ratio 
image (not shown), are generally very flat outside of 
the ionized PAH core. This behavior is consistent 
with the emission expected from VSGs, but is incon¬ 
sistent with the emission properties expected from 
“classical” grains. It also suggests that the PAHs 
and VSGs are well mixed. 



Offset from center (arcseconds) 



Offset from center (arcseconds) 

Fig. 8.— Plots of the observed continuum intensity 
profile of the disk of WL16 along its major axis com¬ 
pared with models of a disk composed only of classi¬ 
cal interstellar grains. Panel a) shows the data and 
model for 10.3 pm; panel b) for 20.8 pm. The solid 
curves represent the observed profiles; the dashed 
curves show the interstellar-grain-only disk model, 
and the dotted curves are the observed PSFs. While 
the disk models are clearly more extended than the 
observed PSFs, they underpredict the observed flux 
from WL16’s disk by an order of magnitude or more 
at radii greater than an arcsecond, and their steep 
slope is incompatible with the broad observed pro¬ 
files. The presence of VSGs is required to account 
for the large spatial extent of non-PAH emission. 

4-3.2. PAH Ionization, Hydrogenation, and Size 
Effects in the WL16 Disk 

The color ratio images presented in Figure 4 can 
give some qualitative indication as to the charge 
state, amount of hydrogenation, and size of the emit¬ 
ting PAH molecules. Quantum mechanical modeling 
and laboratory studies of PAH emission have shown 
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that the charge state of the PAH has an important 
effect on the emitted IR spectrum. In particular, the 
intensities of the C—C stretching and C—H in-plane 
bending modes, which fall in the 6-9 ^m range, are 
generally an order of magnitude stronger in PAH 
cations (singly ionized PAHs) than in neutral PAHs 
(Joblin et al. 1996; Langhoff 1996; Hudgins & Al- 
lamandola 1995a,b). Viewed in this context, it is 
therefore not surprising that the 7.9/8.8 /jm flux ra¬ 
tio shown in the top panels of Figure 4 is essentially 
constant over the entire emitting region of the WL16 
disk, since both of these features (at 7.7 and 8.6 /tin) 
originate from positively charged PAHs. 

DeVito & Hayward (1998) found a central peak¬ 
ing of the continuum-subtracted 8.6/11.3 /irn feature 
ratios in the disk of WL16. We confirm this in the 
central, left-hand panel of Figure 4, which shows 
the 8.8/11.7 ^m flux ratio image. In contrast to the 
featurelessness of the 7.9/8.8 /tm ratio image, a dis¬ 
tinct enhancement within the central ±1.5" core is 
observed (see also Figure 3). 

Recent quantum mechanical calculations show 
that the 8.6 /tm feature originates mostly from C—H 
in-plane bending modes of positively charged PAHs, 
whereas the 11.3 /tm feature originates from C—H 
out-of-plane bending modes of predominantly neu¬ 
tral and anionic PAHs (Bakes et al. 2001). The ob¬ 
served contrast in the 8.8/11.7 /im flux ratio im¬ 
age and cross-cut of Figure 4 may be interpreted 
as a change in charge state of the PAHs from posi¬ 
tively ionized in the central regions to neutral in the 
outer regions (Joblin et al. 1996; DeVito & Hayward 
1998). 

The physical conditions required for the 11.3 /im 
feature to originate mainly from neutral, and not 
positively charged, PAHs, are expressed by the con¬ 
straint, Go/n e < 10 3 , where Go is the UV radiation 
field (6-13.6 eV) in units of 1.6 x 1CV 3 ergs cm~ 2 s _1 
and n e is the electron density in cm -3 (Bakes et 
al. 2001). Of course, the central Herbig Ae star 
is the photo-ionizing and UV source for the WL16 
disk. The calculated UV radiation field in the WL16 
disk at different distances from the central source 
is: G 0 = 1.2 x 10 5 at 0"-l", 2 x 10 4 at l"-2", 
4.8 x 10 3 at 2"-4", and 1.6 x 10 3 for distances at 
4"-6" (DeVito & Hayward 1998). Combining these 
values of the UV flux with the 11.3 /mi feature neu¬ 
trality constraint, leads to derived electron densities 
in the WL16 disk decreasing from n e ~ 100 cm -3 
in the central regions to n e ~ 10 cm -3 in the outer 
regions. 

The bottom panels of Figure 4 show the 12.5/- 
11.7 /in i flux ratios. The most interesting feature 


in this image is the gradual enhancement of the 
12.5/11.7 /tm ratio towards the disk edges. Two 
possible explanations for this enhancement are in¬ 
creasing molecular size and increasing hydrogena¬ 
tion towards the disk edges. 

It has been shown that partially hydrogenated 
small (< 25 C atoms) PAHs, as well as large PAHs 
with 1, 2, or 3 adjacent H atoms can contribute to 
the 11.3 /im feature (Cohen et al. 1985; Mortera & 
Low 1983). However, only relatively large PAHs, 
with 1, 2, or 3 adjacent H atoms, are needed to 
produce the 12.7 /im feature (Duley & Williams 
1981). Quantum mechanical calculations show that 
the 12.7 /im PAH feature strength increases with 
molecule size (Bakes et al. 2001). Therefore, one 
possible interpretation for the enhancement of the 
12.5/11.7 /tm flux ratio at the disk’s edges is the 
growth of larger PAHs (> 50-80 C atoms) at greater 
distances from the central ionizing source. 

Alternatively, since the 11.3 /im feature origi¬ 
nates from PAHs with isolated H atoms (Cohen 
et al. 1985), whereas the 12.7 /tm feature requires 
two adjacent hydrogen atoms on the aromatic ring 
(Verstraete et al. 1996), the slight increase in the 
12.5/11.7 /im ratio apparent towards the edges of 
the WL16 disk might be due to increasing hydro¬ 
genation. Increasing hydrogenation may take place 
due to lower collision rates with electrons that would 
knock away H atoms due to lower n e at the disk 
edges and/or lower UV flux to dissociate H atoms 
from PAHs. 

4-3.3. Disk Asymmetry 

The intensity profiles along the disk’s major axis 
presented in Figure 3 show the presence of an asym¬ 
metry in emission strength with respect to mirror 
reflection about the central position at all observed 
wavelengths. This asymmetry occurs at angular sep¬ 
arations between l"-2"5 from the central position, 
corresponding to radii of 125 AU < r < 300 AU. 
Since this asymmetry is observed at all wavelengths, 
it is likely due to either a column density enhance¬ 
ment in one side of the disk or a gap in the other. 
We do not have enough information to distinguish 
between these two possibilities. 

Although it is tempting to invoke the presence 
of embedded planets to account for the observed 
disk asymmetry due to gravitational perturbations, 
the large radii (125 AU < r < 300 AU) at which 
this asymmetry is observed conflicts with all cur¬ 
rent theories which predict ~ 20 AU as a typical 
outer boundary for planet formation from an accre- 
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tion disk. Furthermore, any planetary time forma¬ 
tion timescale at such large radii would exceed the 
inferred age of the central star. Nevertheless, ob¬ 
jects a few meters to hundreds of meters in size could 
have had time to form at these radii (Kenyon & Luu 
1999). 

4-3.4- Disk Lifetime 

One intriguing problem is the question of how the 
PAH/VSG disk we observe is replenished, since even 
at the disk’s extremities, the timescales for blow¬ 
ing small particles out of the system by radiation 
pressure are only years. Condensation from the gas 
phase is highly unlikely, given the low gas densi¬ 
ties throughout most of the disk, and the high ac¬ 
tivation energies required for PAH synthesis, which 
leads most researchers in the field to conjecture that 
PAHs are formed on the surface of grains. This leads 
us to the suggestion that the most likely scenario 
for replenishment of the PAHs/VSGs in the disk of 
WL16 is from the breakup of larger bodies, such as 
comets or Kuiper belt objects, or from UV or shock 
processing of larger grains as has recently been sug¬ 
gested by Duley (2000). The presence of PAHs in 
a possible planet-forming system has interesting im¬ 
plications for the potential seeding of planetary sys¬ 
tems with organic materials, and is also consonant 
with the discovery of PAHs in meteorites. 

5. Conclusions 

We have presented diffraction-limited images and 
photometry of the disk surrounding the Herbig Ae 
star, WL16, at nine mid-infrared wavelengths, as 
well as photometry of the central star at five addi¬ 
tional near-infrared wavelengths. From our newly 
derived spectral energy distribution for this source, 
we find that the central star and its disk are ob¬ 
served through a foreground screen which provides 
an extinction of Ay = 31 ± 1 magnitudes. We find 
that the star itself has a bolometric luminosity of 
250 L 0 (at an assumed 125 pc distance), and the 
disk contributes « 1% due to very small grain (VSG) 
and polycyclic aromatic hydrocarbon (PAH) feature 
emission. 

The unprecedented spatial resolution of the mid- 
infrared images presented here allows us to confirm 
the presence of ionized PAHs in the central disk 
regions, and we report the discovery of a popu¬ 
lation of larger (> 50-80 C atoms) PAHs and/or 
more hydrogenated PAHs at the disk’s periphery. 
We confirm that a population of graphitic VSGs 
is required to account for the observed large disk 


extent through the five mid-infrared filters which 
exclude the wavelengths of the PAH emission fea¬ 
tures. These VSGs are also required to explain the 
observed mid-infrared excess emission within the un¬ 
resolved inner core (< 0.3") of the WL16 images. 

The disk size is found to be 7"x3.5", correspond¬ 
ing to a cle-projected disk diameter of 900 AU at 
the source. The disk inclination is found to be 
62° ± 2° to our line of sight at P.A. 60°. We find an 
asymmetry in the disk, at all observed mid-IR wave¬ 
lengths, at l"-2"5 from the central source. Finally, 
continuous replenishment of disk material, possibly 
from collisions of larger parent bodies, is required 
to maintain the PAH/VSG disk, which would other¬ 
wise be destroyed by blowout by radiation pressure 
on timescales of years. 
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